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Abstract 

The spectrum of coherent modes for dilute atomic Bose gas, confined in a cylin- 
drical trap, is calculated by applying optimized perturbation theory and the tech- 
nique of self-similar root approximants. The latter technique makes it possible to 
derive accurate analytical formulas. The obtained expressions are valid for arbi- 
trary energy levels of all excited nonlinear coherent modes and for different traps, 
spherical, sigar-shape, and disk-shape. 
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1 Introduction 



At low temperatures, Bose-Einstein condensate of trapped atomic gases is well described 
by the Gross-Pitaevskii equation (see reviews [1-3]). In equilibrium, the condensate state 
corresponds to the ground-state solution of this equation. The mathematical structure 
of the Gross-Pitaevskii equation is identical to the nonlinear Schrodinger equation, which 
in the presence of a confining potential, should possess a discrete set of stationary states. 
Such states can be interpreted [4] as nonlinear coherent modes of trapped atoms [4] cor- 
responding to nonground- state Bose condensates. The excitation of these modes can be 
accomplished by means of resonant modulation of an external potential [4-6]. Such modes 
have also been studied in Refs. [7-9] and a dipole mode has been observed experimentally 
[10]- 

One should not confuse the nonlinear coherent modes, that are stationary solutions to 
the nonlinear Gross-Pitaevskii equation, with elementary collective excitations, which are 
solutions of the linear Bogolubov-De Gennes equations. The nonlinear coherent modes 
are sometimes called topological modes to stress that their spatial behaviour qualitatively 
differs them from each other. Each nonlinear mode generates its own collective excitations 
as small deviations around the given mode. Usually, one considers collective excitations 
above the ground-state mode. These excitations correspond to small density fluctuations 
and are in a reasonable agreement with the Bogolubov spectrum for both superfluid 
helium [11] as well as for trapped Bose-Einstein condensate [12]. In the first case, it is the 
phonon-roton spectrum, which is a unified branch [13,14]. In the case of dilute trapped 
atoms, it is a phonon-single-particle branch. Collective excitations around nonground 
state coherent modes have not yet been considered. 

Before considering perturbations of nonlinear coherent modes, it is necessary to have 
an accurate description of their own properties. In particular, one needs to better un- 
derstand the features of the spectrum of the nonlinear modes themselves. The study 
of this spectrum was initiated in Ref. [4] and, for a weak interaction, was also consid- 
ered in Ref. [9]. In the present paper, we aim at giving a general description of the 
nonlinear-mode spectrum for trapped atoms in a cylindrical trap. We shall describe a 
method for calculating arbitrary energy levels of this spectrum for any magnitude of the 
interaction parameter. We shall also construct accurate analytical expressions for the 
nonlinear-mode spectrum. Having in hands such analytical expressions is convenient for 
studying the spectrum dependence on quantum numbers and system parameters as well 
as for an easier comparison with experiments. 

2 Equation for Coherent Modes 

The interaction of atoms in dilute gases is characterized by the Fermi contact potential 

$(r) = ASM , A = AttH 2 ^- , 

ra 
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in which a s is the s-wave scattering length and mo, the atomic mass. The trapping 
potential is usually taken in the form 



The stationary coherent field <p(r), normalized to unity, is described [3] by the equation 

H[<p(r)]<p(r) = E<p(r) , (1) 
which is often called the Gross-Pitaevskii equation, with the nonlinear Hamiltonian 

H[<p(r)] = -^ + U(r)+NA\<p(r)\ 2 . (2) 

Because of the confining potential U(r), the eigenproblem (1) possesses a discrete spec- 
trum, with the related eigenfunctions being nonlinear coherent modes [4]. 

In what follows, we shall consider a harmonic potential of cylindrical symmetry, with 
a radial frequency 

u±=u x = Uy (3) 

and the anisotropy parameter 

„ S ^. (4) 
It is convenient to work with dimensionless quantities, measuring the space variables 



r\ + r 2 y Tz 



l± 

in units of the transverse oscillator length 



z = ^ (5) 



The coupling parameter 



L = J—. (6) 

V m LOj_ 



= 4vr^ N (7) 
1 1 



is a dimensionless quantity characterizing atomic interactions. Defining the dimensionless 
Hamiltonian and wave function, respectively, 

^^^ft(r), (8) 

in the cylindrical coordinates, we have 

^ = -^V 2 + i(r 2 + ^V)+^| 2 , (9) 
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where 

„ 2 d 2 Id Id 2 d 2 
V = 1 1 1 . 

Then the eigenproblem (1) transforms to the equation 

Hip n = E n ^ n , (10) 

defining dimensionless coherent modes tp n and their spectrum {E n }; n being a multi-index 
labelling the modes. The latter are assumed to be normalized as 



f 

Jo 



roc p2tt f-\-oo 

I r dr dip dz \ip n (r, (p, z)\ 2 = 1 . 

JO J-oo 



Our aim is to find the spectrum {E n } for all energy levels E n of any quantum index n 
and for arbitrary coupling and anisotropy parameters g and u, respectively. 



3 Optimized Perturbation Theory 

To calculate the spectrum of the eigenproblem (10), we shall employ the optimized per- 
turbation theory [15-19], which has been successfully applied to a number of models in 
quantum mechanics, statistical physics, and quantum field theory [15-26]. Numerous 
references on various applications can be found in surveys [27-29]. It is important to 
stress that the optimized perturbation theory was shown to provide accurate results for 
all energy levels and large coupling parameters. 

For the case considered, the procedure can be as follows. We start with the initial 
Hamiltonian of a harmonic oscillator 

H = - l -V 2 + \(u 2 r 2 + v 2 z 2 ) , (11) 

having the oscillator strengths as two trial parameters, u and v. This Hamiltonian pos- 
sesses the eigenvalues 



A = (2n+\m\ + l)u+lj + Z)v, (12) 



nmj \ ■ i" -i ■ -/ - ■ \j • 2 



with the radial quantum number n = 0,1,2,..., asimuthal number m = 0, ±1, ±2, 
and the axial quantum number j — 0, 1, 2, The corresponding wave functions are 



(0) 2n!«H + M 



(n + \m\)\ 

1/4 



•M ( - - ^ 



exp ( - - r ) x 



where L™(-) is a Laguerre polynomial and Hj(-) is a Hermite polynomial. 
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Using a variant of the standard perturbation theory, such as the Rayleigh-Schrodinger 
theory, we may find a sequence {E k (g,u,v)} of the energies E k (g,u,v) for the approxi- 
mation orders k — 0, 1, 2, Here, for the simplicity of notation, the dependence of the 

energy levels on the quantum numbers is not written down explicitly, but it is assumed. 
For instance, E (g,u,v) = E^j. Then the trial parameters are to be transformed to 
control functions [15-19] u k (g) and Vk(g) such that to render the sequence {e k (g)} of the 
optimized terms 

e k (g) = E k (g, u k (g),v k (g)) (14) 

convergent. The control functions can be found from optimization conditions. For exam- 
ple, we may employ the fixed-point condition 



' d d\ 
Su du +Sv dvl Ek ^ 9, V ' = ° ' 



(15) 



whose solutions are u = u k (g) and v = v k (g). The optimized approximants (14) are 
valid for arbitrary quantum numbers and for the whole range of the coupling parameter 
g. The same procedure can be used for anharmonic traps, described by anharmonic 
confining potentials of different powers [29], integer or noninteger. Note that the usage of 
anharmonic traps may be important for some experiments [30] dealing with elementary 
collective excitations. 

In the first order, we have the energy 



E 1 (g, U ,v) = P -(u + 1 -) + l 



v + 



+ uy/v gL, 



in which the notation 

p = 2n + |m| + 1 , q = 2j + 1 
for the combinations of quantum numbers is introduced and where 



(16) 
(17) 



L nmj 



= ~^1= I l^nljir, <P, z)\ A rdr dy? dz 

U\/V J 



7T^ 



(n + \m\)\ 2ij\ 



lyjv 

f 
Jo 



x 2\m\ e -2x 



X 



dx 



r 

Jo 



The fixed-point condition (15) gives the equations 

p(i-A) + -#=°. ■ 

V u z J pv\ q 




e~ 2t Hf(t) dt . 



= 



(18) 



pv\jvq 

for the control functions u = U\(g) and v = V\(g), where, for convenience, the notation 

s = 2p^/q g vlnmj (19) 

is used. Substituting the solutions for the control functions, given by Eqs. (18), in the 
form (16) results in the optimized approximant e±(g), according to definition (14). This 
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procedure can be continued to a desired approximation order. Here we limit ourselves by 
the first-order approximation 

E = e 1 (g) = E 1 (g,u 1 (g),v 1 (g)), (20) 

in which the control functions Ui(g) and Vi(g) are the solutions of Eqs. (18). The latter 
equations, for a given set {n,m,j} of quantum numbers and a coupling parameter g, 
require a numerical solution. Hence, the spectrum (20) can be calculated only numerically. 
However, it is always desirable to possess an approximate analytical expression that could 
be easy to treat with respect to varying system parameters, such as quantum numbers 
and the coupling parameter. 



4 Self- Similar Root Approximants 

In order to derive analytical formulas for the spectrum of the nonlinear coherent modes, 
we may employ the technique of self-similar crossover approximants [31,32]. For this 
purpose, we need to know the asymptotic expansions of the spectrum in the limits of the 
weak and strong coupling. These expansions for the spectrum (20) can be derived from 
Eqs. (16) and (18). It is convenient to use the variable (19) that is proportional to the 
coupling parameter; so that if g — > 0, then s — ■> 0, and when g — > oo, then s — > oo. In the 
weak-coupling limit, we find the expansion 

E ~ a + ais + a 2 s 2 + a 3 s 3 , (21) 

as s — > 0, where 

qv 1 p + 2qv (p + 2qv) 2 

«o=P + y, 01 = 2 p(quy/ 2 ' ° 2 = ~ 16p 3 (qu) 2 ' ° 3 = Up^qv) 7 / 2 ' 

In the strong-coupling limit, we obtain 

E ~ M 2/5 + M~ 2/5 + hs-^ 5 + b 3 s- 2 + M- 14/5 + M~ 18/5 , (22) 
as s — > oo, where 

Ab = 5 , 46i = 2p 2 + (qv) 2 , 20b 2 = -3p 4 + 2p 2 (qu) 2 - 2(qu) 4 , 

20b 3 = 2p % - p 4 (qu) 2 - 2p 2 (qvf + 2(qvf , 

5006 4 = -44/ + 22p 6 (qu) 2 + 2p 4 (qu) 4 + 78p 2 (quf - 69(qu) 8 , 

125006 5 = H22j9 10 - 595p 8 (qu) 2 - 70p 6 (qu) 4 + 440p 4 (^) 6 - 3640p 2 (qu) 8 + 2821(gz/) 10 . 

The asymptotic expansions (21) and (22), valid in the weak-coupling limit s — > and, 
respectively, in the strong-coupling limit s — > oo, can be sewed by using the technique of 
self-similar crossover approximants [31,32]. This technique makes it possible to construct 
interpolative formulas that provide correct asymptotic expansions, at the same time giving 
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good accuracy in the whole interval of the variable s. The detailed description of the 
method can be found in Refs. [31,32]. Taking into account only one term in the strong- 
coupling expansion (22), we get the first-order self-similar root approximant 



E{ = oo(l + As) 



2/5 



(23) 



in which 



a = p + 



qv 



Aal 12 



1.746928 



Retaining two terms in Eq. (22) yields the second-order approximant 



E* 



a 



[l + A lS f 5 + A 2 s 



where a is the same as in Eq. (23) and 



A ia f /6 = 2 o33913 
Similarly, in the third order we find 

E; = a {[(1 + B lS f 5 + B 2 s 2 
where ao is again as earlier and 



2p 2 + (qvf 



5/6 



1/5 



A 2 al = 3.051758 



(24) 



11/10 



+ B 3 s : 



;| 2/15 



(25) 



R 125/22 



2p 2 + (quf 



5/66 



= 1.405455 

10/11 



8p 4 + \2p 2 (qv) 2 + (qu) 4 



5/6 



B 3 a 5/2 = 5.331202 . 



B 2 a 7 Q /n = 6.619620 [2p 2 + (qu) 2 

In the same way, we may construct higher-order approximants, which, however, we will 
not write down explicitly. 

Let us stress that the self-similar root approximants (23) to (25) are valid for the whole 
range of variable (19), related to the coupling parameter (7); they are also valid for an 
arbitrary anisotropy parameter (4), incorporated in the variable (19), as well as for all 
quantum numbers that enter expressions (23) to (25) through the parameters (17) and 
the variable (19). The accuracy of the root approximants can be checked by comparing 
their values with those obtained numerically from Eqs. (18) and (20). 

Bose-Einstein condensates are often treated by using the Thomas-Fermi approxima- 
tion, when one omits the kinetic energy operator in Eq. (10). In such a case, the solution 
to Eq. (10) reads 



ij TF (r,z) = (r 2 - r 2 - v 2 z 2 ) 



2 2\ I/ 2 



2g 



where Tq = 2E?f and the energy is obtained from the normalization condition for i^tf, 
which results in the expression 



TF 



15 n 2 / 5 
— qu 
4tt J 
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The Thomas-Fermi approximation is not self-consistent since the average of the Hamil- 
tonian (9), calculated with the wave function ipTF, diverges. The energy Etf becomes 
correct for an asymptotically large coupling parameter g — > oo, but for small and inter- 
mediate g this formula is not accurate. There have been several suggestions of improving 
the Thomas- Fermi approximation by means of additional corrections [33-35]. However, 
even after being improved, this approximation provides us solely the ground-state energy, 
being enable to describe the excited energy levels of the nonlinear coherent modes. 

Figure 1 shows the ground-state energy of the condensate in a cigar-shape trap, with 
v = 0.1, calculated in different approximations. For comparison, the Thomas- Fermi 
energy is also presented although its accuracy up to the coupling parameters g ~ 300 is 
quite bad. The maximal errors for the self-similar root approximants from E* to are, 
respectively: 8%, 3.5%, 2%, 1.2%, and 0.8%, which demonstrates good convergence. 

The accuracy of the self-similar root approximants for different energy levels of the 
coherent modes, varying the coupling parameter g, and for different trap shapes is illus- 
trated in Figures 2 to 4, where the percentage errors as functions of g are drown for a 
cigar-shape trap (Fig. 2), spherical trap (Fig. 3), and a disk-shape trap (Fig. 4). In all 
cases, there is the following standard situation. The maximal, with respect to g and 
quantum numbers, percentage errors for E{ are between 4 — 12%; for E^, between 2 — 5%, 
and for £3, of order 1%. The approximants E\ and E\ are close to £3, because of which, 
not to overload the Figures, they are often omitted. 

5 Conclusion 

The nonlinear coherent modes of Bose-Einstein condensate at zero temperature are stud- 
ied. These are described by the Gross-Pitaevskii equation and correspond to the sta- 
tionary solutions of the latter. In the presence of a trapping potential, the spectrum of 
the coherent modes is discrete. Our main aim in this paper was to find a convenient 
and accurate way of calculating the energy levels of this spectrum for arbitrary quantum 
numbers, any values of coupling parameters, and for different aspect ratios of a cylindri- 
cal trap. One way for calculating the spectrum is by employing optimized perturbation 
theory, which requires numerical solution of the equations for control functions. Another 
possibility is to invoke the technique of self-similar root approximants, which results in 
sufficiently simple, and at the same time accurate, analytical formulas for the spectrum. 
The advantage of possessing analytical expressions is in the feasibility of a relatively easy 
study of their behaviour with respect to all parameters characterizing the system. And 
to understand better the properties of the coherent modes is necessary for choosing the 
optimal conditions for their experimental observation and practical usage. 
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Figure Captions 



Fig. 1. The ground-state energy for a cigar-shape trap, with v = 0.1, as a function of 
the coupling parameter g, calculated by using formula (20) (solid line), and the self-similar 
root approximants E{ (long-dashed line), E\ (short-dashed line), and E\ (dotted line). 
The energy Etf in the Thomas-Fermi approximation is shown by the dashed-dotted line. 

Fig. 2. The percentage errors of the self-similar root approximants for the energy 
levels of coherent modes in a cigar-shape trap, with v — 0.1, corresponding to E\ (solid 
line), E% (long-dashed line), E% (short-dashed line), El (dotted line), and E^ (dashed- 
dotted line) for different quantum numbers: (a) n = m = j = 0; (b) n — j — 0, m — 2; 
(c) n — j — 0, m — 10. 

Fig. 3. Percentage errors of the self-similar root approximants for E{ (solid line), E^ 
(long-dashed line), and £g (short-dashed line) for a spherical trap, with v — 1, and for 
different energy levels: (a) n = m = j = 0; (b) n — j — 0, m — 2; (c) n — 3, m — 2, j = 
1. 



Fig. 4. The case of a disk-shape trap with v = 10. Percentage errors of the ap- 
proximants El (solid line), E\ (long-dashed line), E£ (short-dashed line) and E\ (dotted 
line) for different coherent modes: (a) n = m = j = 0; (b) n — j — 0, m — 2; (c) 
n — 3, m — 2, j — 1. 
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